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A large reduction of the p° mass in the nuclear medium is reported, inferred from dipion pho- 
toproduction spectra in the 1 GeV region, for the reaction _ffe(7,7r + 7r~)X with a 10% duty fac- 
tor tagged-photon beam and the TAGX multi-particle spectrometer. The energy range covered 
(800<E 7 <1120 MeV) lies mostly below the free p production threshold, a region which is believed 
sensitive to modifications of light vector-meson properties at nuclear-matter densities. The p° masses 
extracted from the MC fitting of the data, m* = 642T40, 669T32, and 682±56 MeV/c 2 for E 7 in 
the 800-880, 880-960, and 960-1040 MeV regions respectively, are independently corroborated by a 
measured, assumption-free, kinematical observable. This mass shift, far exceeding current mean- 
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field driven theoretical predictions, may be suggestive of p decay within the range of the nucleonic 
field. 

I. INTRODUCTION 

Hadronic dynamics have until recently been comprising of two non-overlapping domains, distinctly separated along 
the lines of their respective description of matter as hadronic or quark. On the low-energy side of the spectrum, 
matter is probed on a scale where pions and baryonic resonances are the relevant constituents. These have been 
rather succesfully employed by quasi-phenomenological models in providing prescriptions for free and in-medium 
hadronic interactions. At the higher-energy end, far above energies typical of baryonic resonances, quark degrees of 
freedom become accessible and matter tends asymptotically towards quark-gluon plasma. In this regime, described 
by perturbative QCD, quarks become deconfined and chiral symmetry is restored. 

Although the origin of quark confinement is not known, it is evident that, whatever the reasons, it must induce 
spontaneous breaking of chiral symmetry The two phenomena are therefore linked and in the limit of chiral 
symmetry restoration, as quark masses tend to zero, vector- meson masses and widths are also expected to change [p|. 
The transition from hadronic to quark matter and the effect of increasing density and temperature on the properties 
of light vector mesons have been addressed in the context of QCD sum rules jjj as well as effective Lagrangians 
0. In particular, a temperature and density dependent lowering of the p° mass is regarded as a precursor of the 
chiral phase transition, expected to be measurable even at normal nuclear density Q . Thus, the energy region from 
hadronic phenomenology to the domain of perturbative QCD has come to the foreground, as the study of vector- 
meson property modifications in this range appears to be holding the key to understanding the mechanism of chiral 
symmetry restoration. 

The recent advent of high duty-cycle photon beams, on one hand, and high-resolution dilepton spectrometers, on 
the other, have opened up the possibility of reconciliation between the hadronic and quark depictions of matter. 
Hadronic probes coupled with hadron spectrometers had been extensively used in the past for energies up to the 
A resonance, but higher-up, where multi-pion production channels increasingly dominate, the combination becomes 
rapidly cumbersome due to medium distortions from initial and final state interactions. For the study of the vector 
mesons that couple to multi-pion states, the photon is ideally suited. The recognition of the electromagnetic interaction 
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as an indispensable probe at either the entrance or the exit channel, if the vector-meson field in the nuclear medium 
is to be delineated, has generated activity in a variety of fields. 



A. Relativistic Heavy Ion Results 

At relativistic energies, under extreme conditions of temperature and density, the quest for signs of a phase transition 
from hadronic matter to quark-gluon plasma is at the heart of experimental programs using heavy-ion beams (e.g. 
at CERN, GSI, and RHIC) and dilepton spectrometers. A series of pioneering studies from the CERES, HELIOS-3, 
and NA50 colaborations at SPS/CERN with central S + Au, S + W, Pb + Au, and Pb + Pb collisions, have largely 
been interpreted as indicative of a downward shift of the p° mass in the nuclear medium [|]j7j . The invariant mass 
spectra for dilepton production that have been measured in these experiments, when compared with the respective 
p-A spectra, show a large enhancement at low invariant mass regions. The excess dileptons are thought to originate 
from the decay of mesonic resonances produced in ir + ir~ annihilation and, in the mass region 0.2<m;+/- <0.5 GeV/c 2 , 
the models advocating vector-meson medium modifications attribute the enhancement to the decay of the p° meson 
with a downward-shifted mass Nonetheless, an explanation in terms of conventional phenomenological p°-meson 
medium modifications is also consistent with the CERN data |l^,|ll|,|l2| . In this more conservative approach, the 
p spectral function below 0.6 GeV/c 2 in dilepton invariant mass is appreciably enhanced from the contributions of 
"rhosobar"-type excitations such as the A7V -1 , N*(1720)N^ 1 and A*(1905)-/V _1 , a consequence of the strong coupling 
of the p° meson with it + tt~ states in the nuclear medium. This being the case, the downward shift of the p -meson 
mass in the nuclear medium may amount to no more than a convenient parametrization 

From a theoretical perspective, a model combining chiral SU(3) dynamics with vector-meson dominance in an 
effective Lagrangian has shown that chiral restoration does not demand a drastic reduction of vector-meson masses 
in the nuclear medium p"3| ] . The latter model, in qualitative agreement with the hadronic-phenomenological models 
p0| , p2| , predicts a substantial enhancement of the p° spectral density below the nominal resonance mass, with only 
a marginal mass reduction. This result is tantamount to the p° dissolving in the medium. Moreover, both an 
interpretation of the CERES data as the outcome of either medium-modified hadronic interactions, or interactions of 
dissociated quarks in the quark-gluon phase, yield remarkably consistent results in the framework of Ref. [ft3f , leading 
to the conjecture that both the hadronic and quark-gluon phases must be present |fL3| . The latter conclusion is also 
drawn by a synthesis of the mass-scaling Q and hadronic-phenomenological [jl2| models, leading to the prediction 
that dynamical-hadronic effects that are dominant up to about nuclear density, mainly via the highly-collective 



N*(1520)N _1 state, gradually give way to p°-meson mass scaling as the quark degrees of freedom become increasingly 
relevant. In the intermediate crossover region of "duality", both the hadronic and quark-gluon phases of matter are 
expected to coexist O]. Thus, although a consensus has been reached on the issue of coexistence of the hadronic 
and quark-gluon phases in the transition region roughly placed upwards from ^1 GeV, the question remains as to the 
extent to which the nuclear medium induces modifications to the p° mass and width. 

Other less directly related experimental conjectures of medium modifications of the p° mass have been deduced from 
an anomalous J/V> suppression, reported by the NA50 collaboration for Pb + Pb collisions Jl5| ], enhanced K + - 12 C 
scattering cross sections fltst] , and an IUCF experiment of polarized proton scattering on 28 Si, though in the latter the 
medium renormalization of the p° mass required for agreement with the data is inconsistent for different observables 
measured in the same experiment [fl7f . 

In summary, the experimental results discussed so far are inconclusive with regard to the magnitude of the p°- 
meson medium mass modification, and underline the limitations encountered in complex processes, where the probe 
interacts with nuclear matter via the strong interaction and the channel being investigated may not be disentangled 
from conventional medium effects. 



B. Electromagnetic Probe Results 

These difficulties are largely overcome with the use of photons which do not suffer from initial state interactions. 
The availability of high-flux photon beams, which have compensated for the low interaction cross sections with 
nuclear matter, complemented by wide-angle multi-particle spectrometers, have made possible a new generation of 
experiments. In the E 7 ~1 GeV region, matter is probed at short distances <1 fm, which is at the gateway of the 
energy range where vector-meson properties are expected to undergo modifications. In this domain, baryons and 
mesons are still the relevant constituents for the description of matter, yet their quark content becomes increasingly 
manifest, a fact which is reflected in QCD-inspired phenomenological models |l8|]l9| ]. The p° meson is the best 
candidate among the light vector mesons p 0,± , ui, and ip as a probe of medium modifications, since, due to its short 
lifetime and decay length (1.3 fm), a large portion of p° mesons produced on nuclei will decay in the medium. 

For effective measurements with the photon as probe in the 1 GeV region, kinematically complete experiments 
are required, which in turn necessitate the use of high-duty-cycle tagged photon beams and large-acceptance multi- 



particle spectrometers. These requirements were met by the TAGX detector |20 at the Institute for Nuclear Study 
(INS), where the TAGX collaboration has completed a series of experiments with the 3 ffe(7,7r + 7r~)X reaction and 



a 10% duty-cycle tagged photon beam. First, a lower-energy experiment (380<E 7 <700 MeV) measured the single- 
and double- A contributions to tt + it~ production Having established these important non-p dipion processes, 
the kinematics of the 3 i7e(7,7r + 7r~)X reaction were investigated in the range 800<E 7 <1120 MeV, aiming at the 
p° — > tt + tt~ channel. 

Both the coherent and quasifree p° photoproduction mechanisms are relevant in the energy region of interest. For 
a 3 He target, the energy threshold for the former is E 7 m 873 MeV, whereas the 1.083 GeV energy threshold of the 
elementary p° photoproduction reaction on a nucleon is lowered for the quasifree channel in the nuclear medium, as 
the Fermi momentum of the struck nucleon may be utilized to bring the p° meson on shell. Coherent photoproduction 
on nuclei is characterized by small four-momentum transfers, resulting in p° mesons which mostly decay outside of 
the nucleus |22|j . The latter is, consequently, of limited utility in probing vector-meson medium modifications. 

In contrast, quasifree subthreshold photoproduction (E 7 < 1.083 GeV) on one hand warrants that the interac- 
tion took place in the interior of the nucleus, since the nucleon Fermi momentum is required to produce the p°, 
and on the other produces slower p° mesons, more likely to also decay inside the nucleus. Moreover, subthreshold 
photoproduction, whether the target nucleus remains bound in the final state (exclusive channel) or decomposes to 
its constituent nucleons (breakup channel), may be correlated with either a nuclear mean- field or nucleonic medium 
effect on vector- meson properties (Section |VB| ). Specifically, exclusive subthreshold photoproduction produces p° 
mesons which, on the average, traverse distances comparable to the size of the nucleus before their hadronic decay. 
This process is therefore a probe of vector-meson medium modifications at normal nuclear densities, a regime which 
has been the focus of mean-field driven theoretical models. The breakup channel, on the other hand, is more likely 
to produce p 0, s which are slower relative to the struck nucleon and may travel distances shorteer than the nucleonic 
radius before decaying (Section |VB| ). Thus, large densities in the interior of the nucleon may become accessible via 
the breakup channel, amounting to a nucleonic medium effect. This is the realm of the emerging hadronic-quark 
nature of matter, a domain virtually unexplored. In the deep subthreshold region, and for large momentum transfers 
to the target nucleus, the coherent and exclusive-quasifree channels are suppressed, and the breakup-quasifree process 
becomes dominant. It is the subthreshold dynamics that motivated the 800<E 7 <1120 MeV p° photoproduction 
experiment. 

The aim of the experiment being as stated above, the choice of 3 He as the target becomes almost ideal. The low 
photon energies utilized to induce subthreshold photoproduction result in slow p° mesons with a small Lorentz boost, 
and therefore a large probability for decay within the nuclear volume. Thus, in the case of the exclusive channel, 
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a large nuclear radius is not necessary, and the larger nuclear density of a heavier target is predicted to have only 
a marginally enhanced effect on the p° meson mass Jf9[ . Furthermore, if the breakup process is dominated by the 
nucleonic field, the size of the target nucleus is irrelevant. Finally, the 3 He target is the lightest nucleus where a 
nuclear medium effect may be discernible, without the complexity of overwhelming final state interactions (FSI). 

The p° detection is further aided by the inherent selectivity of the TAGX spectrometer (Section [II B| ) to coplanar 
7r + 7r~ processes |20j. This is due to the limited out-of-plane acceptance of the spectrometer, which preferentially 
selects the p° — > ir + ir~ channel, at the expense of two-step processes (e.g. FSI) and uncorrelated ir + ir~ production 



at distinct reaction vertices, the latter accounting for the majority of non-/?° background events (Section IV A ). This 
favorable feature of the spectrometer promotes an otherwise small component of the total amplitude, namely the 
subthreshold breakup channel, to a sizeable experimental signal. 

The mass of the p° meson in the nuclear medium was extracted from the dipion spectra of the 800<E 7 <1120 
MeV experiment with the aid of Monte Carlo (MC) simulations p3) . The reported mass shift was far larger than the 
predictions of any mean-field driven model pertaining to the exclusive process for 3 He [ fl9| ; a calculation based on QHD 
assuming a deep scalar potential yielding a p°- 3 He bound state, on the other hand, produced much better agreement 
p4j . The result of Ref. |2j| led to a reanalysis of the lower-energy measurements [^lj including the p° — > tt + tt~ 
channel, and allowing for p° production with a shifted mass. The outcome of this reanalysis was an even larger shift 
p5fl , possibly indicating a mechanism other than those previously considered. A nucleonic medium effect, as sketched 
earlier, may be consistent with a large p°-mass reduction, although a theoretical model has yet to be fully developed 
p6jp7| . Though more work is needed to firmly establish and better understand these results pgf] , the 3 i/e(7,7r + 7r~)X 
experiment for photon energies in the range 800<E 7 < 1 120 MeV constitutes the first direct measurement of the p° 
mass in the nuclear medium. In this report, a new and more thorough analysis is presented, including the first direct 
evidence of the characteristic J = 1 signature of the p° meson decay in the subthreshold region, as well as refinements 
in the simulations and fitting procedure, relative to the analysis of Refs. |2T|]23| , |2l| , leading to higher confidence in 
the extraction of the in-medium p° mass. 

The paper is organized in six sections. In Section [n], the experimental set-up and the calibration procedure are 



reviewed. In Section [II the data-analysis algorithm is outlined in conjunction with the experimental aspects of 
Section y. Section [jv| focuses on the MC techniques and the fitting of the data, leading to the extraction of the mass 
shift. In Section [v| the data are compared with the MC calculations, and the p° mass shifts are discussed. Finally, 



the conclusions are presented in Section VI. 
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II. APPARATUS 



The INS tagged photon beam and the TAGX spectrometer (Fig. |l|) ^0|, the new straw drift chamber |2(j[j , and 
different aspects of the data-analysis procedure have all been described in detail elsewhere, where the reader 

is referred for a more extensive discussion. In this section, a brief overview of the experimental apparatus is provided, 
as it applies for the 3 £/e(7,7r + 7r~)X measurements, stressing the elements that were either introduced for the first 
time in this experiment, or that are important for the data analysis. 

A. Photon beam and A He target 

The photon beam is produced utilizing the 1.3 GeV Tokyo Electron Synchrotron. A scries of innovative technical 
improvements led in 1987 to the upgrading of the photon beam to one of medium duty cycle. In the present experi- 
ment, the endpoint electron energy E s is 1.22 GeV at an average 10% duty factor, corresponding to a 5 ms extraction 
time. The instantaneous energy of the extracted electrons has a sinusoidal dependence, and it is known by measur- 
ing the extraction time. The extracted electrons are directed via a beamline onto a thin platinum radiator where 
bremsstrahlung photons are produced, while the scattered electrons are bent away from the beam by a rectangular 
analyzer magnet of 1.17 T. The magnet settings of the extraction beamline vary sinusoidally in time, in phase with 
the E s energy. An array of 32 scintillator electron-tagging counters, each with a 10 MeV/c momentum acceptance, 
detect the scattered electrons. The position of the tagger registering the scattered electron determines its energy, and 
consequently that of the bremsstrahlung photon as E 7 = E e - E e < with AE 7 ~ ±5 MeV accuracy. A second set of 
8 backing taggers participates in the coincidence triggering signal, along with the 32 frontal ones, and is discussed 
in the following section. The tagged photon intensity was maintained at an average of ~3.5 x 10 5 7/s, well within 
the tolerance of the data acquisition system for accidental triggers. The photons, distributed over a beam spot of 
^2 cm in diameter, are subsequently incident on a liquid 3 He target. The target is at a temperature of 1.986±0.001 
K, corresponding to a 0.0786 g/cm 3 density, and is contained in a cylindrical vessel 90 mm in height and 50 mm in 
diameter at the center of the TAGX spectrometer [ ^0[ . 

The tagger hits are related to the photon flux incident upon the target after efficiency corrections. In particular, 
due to the collimation of the photon beam downstream from the taggers, some of the tagged photons do not reach 
the target. To determine the tagging efficiency, and consequently the photon flux, a lead-glass Cerenkov counter is 
placed in the photon beam, with reduced flux, downstream from the target and the tagger scalers are periodically 
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recorded both with and without the platinum radiator in place. The efficiency per tagger counter is determined by 
the relation 
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where T l is the scaler count for each of the frontal taggers. The term 
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corresponds to the coincidences 



between a tagger-counter and a lead-glass Cerenkov hit from a spurious photon not originating from the platinum 



radiator, and 
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is the number of accidental coincidences with the radiator in place, but with the Cerenkov 



hit registering with a delay of the order of 100 ns with respect to the tagger counter signal. These two terms turn 



out to be negligible relative to the term 
photons per tagger reaching the target as 
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The efficiencies r] t and radiator out/in ratios TZ t — T iout /T Hn for the 32 frontal taggers are recorded in a number 
of dedicated runs, in regular intervals throughout the experiment, and they are used, along with Eq. (^), in the 
empty-target background subtraction procedure applied to the measured spectra (see Section III). 



B. Spectrometer 

The TAGX spectrometer has an acceptance of 7r sr for charged particles (neutral-particle detection was not utilized 
in the present experiment), and has been in use at INS since 1987. It consists of several layers of detector elements 
(Fig. |l]) positioned radially outwards from the target vessel, which is located at the center of the 0.5 T field of a dipolc 
analyzer magnet. 

Directly surrounding the target container is the inner hodoscope (IH) , made up of two sets of six scintillator counters, 
one on each side of the beam. The IH is used in the trigger, as well as in measuring the time of flight (TOF) of the 
outgoing particles (2^j2^]. As it is placed inside a strong magnetic field, the light signal is carried by optical fibers to 
the photomultiplier tubes, which are located at the fringes of the magnetic field two meters away. 

Next is the straw drift chamber (SDC), operating since 1994 and installed expressly for the measurement of the p° 
mass in 3 He, with the objective of improving the momentum resolution for the detection of the p° — ► it + tt~ decay 
channel p9[ ]. The SDC consists of two semi-circular cylindrical sections, each containing four layers of vertical cells. 
The "straw" cells have tube-shaped cathodes which induce a radial electric field, and consequently have a regular field 
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definition and high position resolution (~150 /im). The SDC was designed to preserve the 7r-sr large acceptance prior 
to its installation, to not impose extensive modifications of the spectrometer, and to not induce significant energy 
losses to traversing particles by keeping its thickness to minimum. The installation of the SDC required, nonetheless, 
the replacement of the IH from an earlier set of scintillators with the one described above. 

Surrounding the SDC are two semicircular cylindrical drift chambers (CDC) subtending angles from 15° to 165° on 
both sides of the beamline in the horizontal plane, and ±18.3° in vertical out-of-plane angles. The CDC is composed 
of twelve concentric layers of drift cells, yielding a ^250-300 /im horizontal and 1.5 cm vertical resolution. Together 
with the SDC, they are used to determine the planar momentum and emission angle of the charged particles traversing 
them, and the vertex position of trajectory crossings. 

The outer set of 33 scintillator elements comes next, serving as the outer hodoscope (OH), with PMT's attached at 
both ends to help determine the track angle relative to the median plane. The two sets of hodoscopes, IH and OH, 
measure the TOF corresponding to the tracked trajectories. 

Other components of the TAGX spectrometer include 4 sets of 155 mm x 50 mm x 5 mm e + e~ background veto 
counters positioned along the OH arms in the median plane. The veto counters eliminate charged-particle tracks 
registering within Az = ± 2.5 mm, mostly affecting forward- focused e+e - pairs produced copiously downstream from 
the target, but having a small effect on tt + tt~ events. 

C. Data acquisition and calibration 

The channel of interest being tt + tt~ production from the decay of the p° meson, two-charged particle coincidences 
on opposite sides of the beam axis were required of the trigger. Two levels of triggering are implemented in order to 
optimize the data acquisition electronics. The pretrigger 
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VT = IH L -IHr-Y, TAG lback ■ EM front (3) 
»=i 

is generated within 100 ns from the occurence of an event. A coincidence of a left and right IH hit with a backing 
tagger hit is required, and not rejected by the forward e + e~ veto counters. The main trigger 

32 

MT = VT- OH L ■ OH R ■ TAG lfront ■ hMMt (4) 

i=i 

requires the coincidence of the pretrigger with a left and right OH hit and a forward tagger hit, not rejected by the 
computer inhibit signal. A window of 400 ns is available between the VT and the MT, after which the CAMAC is 
cleared for the next VT. Typical counting rates are 2 kHz and 30 Hz for the VT and MT, respectively. 
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The calibration of the scintillation counters and the CDC and SDC have been the subject of extensive effort 
p5] , p9|j3l| . More recently, a series of modifications implemented in the track fitting algorithms has resulted in significant 
improvements, mainly in the planar-momentum resolution |32]. It is the tracking that is discussed next. 

The CDC consists of four groups of three- wire layers (Fig. [j]). The last layer of wires for each group was intended 
for charge division readout and had not been employed in the past p0| . Instead, hits from the eight remaining CDC 
layers were used for the reconstruction of the planar momentum p xy , emission angle 9 SC , and vertex position (see 



also Section III B| ) . This earlier eight-layer tracking procedure did not incorporate the SDC information either, thus 
resulting altogether in a less-than-optimal momentum resolution. Since longer effective lengths of reconstructed tracks 
result in higher-quality fits, however, TDCs from the last layer of wires of the fourth group has been implemented 
for the first time in the present analysis. Furthermore, the SDC data have also been used for the first time, a 
combination which yields an overall improvement in the planar momentum resolution estimated to be <J Pj:y /p X y — 
0.0892pj;j, + 0.0057, compared with o Vxy j\>xy — 0.1150pj;y + 0.0078 from the 8-layer CDC analysis [32| . For a particle 
of p £C j / =300 MeV/c, this amounts to a 40% improvement in the planar momentum resolution. The corresponding 
improvement in the vertex position is reflected in Fig. ^, and has allowed for more stringent tests in the selection of 
two-track events which originate from the target area. A minor improvement has also been noted in the emission 
angle resolution, which stays relatively constant at a v ~ 0.3° throughout the range of typical planar momenta 100 
MeV/c < p xy < 500 MeV/c @. 

The steps involved in the tracking of trajectories through the SDC and CDC may be summarized in the following: 
The CDC TDCs operate in "common-stop" mode, with the start determined from each CDC sense wire and the stop 
from the IH ^cj. The CDC drift times are first corrected by the corresponding TDC timing offsets. The drift-length 
to drift-time relation is determined next, per layer of CDC wires, as a fifth-order polynomial. This is an iterative 
process, where an initial set of parameters is used to reconstruct a sample of well-defined high-momentum tracks. 
The reconstructed trajectories yield a new set of parameters, and the procedure is repeated until the convergence 
condition is reached, namely that the residual root mean square (RRMS) improvement over the final two cycles is no 
better than 0.5%. 

Once a CDC track has been reconstructed, a similar procedure is followed for the SDC, where first the TDC timing 
offset is determined, and subsequently a SDC length-to-time relationship is extracted. This accomplished, "best" 
SDC tracks are identified, which qualify as candidate extensions of a selected CDC track. Typically, 2-4 SDC tracks 
are selected as possible extensions of a reconstructed CDC track when all four SDC layers have registered a hit, to 
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a maximum of 8 candidate tracks if one SDC layer is missed. The SDC tracks are approximated by straight-line 
segments, since the error in the position of even the slowest particles which may be expected to result in valid two- 
track events is within the 150 /im tolerance of the SDC. Finally, the SDC candidates are matched with the CDC 
reconstructed track, by requiring the minimal CDC + SDC RRMS of the combined track. 

The obtained TOF resolution a t is better than 380 ps pOfl . The TOF is used for particle identification, as well as 
for the determination of the particles' OH position (along the z-axis in the TAGX frame as shown in Fig. Ill) . 



III. DATA ANALYSIS 

The data presented in this report were collected in two periods, with 3 He and empty-target measurements in each. 
The superior quality of the photon beam, and a longer running period, resulted in better statistics and a higher 
ratio of ir + n~ to accidental triggers for the second phase. This is reflected in the tagger efficiencies and radiator 
out/in ratios, defined in Eq. (Q), which during the later part of the experiment were generally improved. A total of 
16,366 7r + 7r~ events have been identified from the analysis of two-track events, comprising 73% of the total number 
of reconstructed events, the remaining being of three (23%) or more tracks (<4%). With the extraction of the ir + ir~ 
yield Y, the total cross section is determined from the relation 

Y 

<jt = ,, (5) 

N T N 7 T] W + W -T] daq 

where Nt (nuclei/cm 2 ) is the 3 He target density seen by the photon beam, and iV 7 is the incident photon flux. 
The efficiencies rjdaq and 77^+^- account for the data-acquisition livctime and 7r + 7r~ detection efficiencies. The latter, 
which is in the range of 2.7-6.8% for the p° channels, is determined by dedicated MC routines and is reaction-channel 
specific (see Ref. |20|). 



A. Empty target background 



In Section |II At the extraction of the tagger efficiencies, and the normalization of the tagger scalers to reflect the 
number of photons incident on the target, were discussed. These are utilized in determining the appropriate factor 
by which empty-target spectra are scaled prior to their subtraction from the corresponding 3 7Je-target spectra, for 
the removal of target background counts. The procedure is briefly summarized in the following steps. 

At regular intervals throughout each of the 3 ffe-target and empty-target running periods, the lead-glass Cercnkov 



counter is employed in dedicated reference runs to determine the quantities ^ and as described in Section [I A 
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The total number of photons incident on the target per experiment is extracted as the sum of the raw scaler counts 
T,, recorded for each run, corrected by the efficiency and out /in ratios for that run, according to Eq. (0), and weighted 
by the normalized energy distribution of the scattered electrons. In particular, 77, and lZ t for each run are calculated 
from the corresponding quantities of the reference runs, on the assumption that they vary linearly with the raw 
scaler counts accumulated between runs. The ratios of photon fluxes between each 3 f/e-target and its corresponding 
empty-target experiment yield the scaling factors by which the latter are normalized prior to subtraction from the 
former. The x-coordinate spectrum of the two-pion crossing vertex after background subtraction, indicative of the 
accuracy of this procedure, is shown in Fig. |^. 



B. Experimental observables 



The calibration procedure, discussed in Section II C, allows the extraction of the planar momentum p xy , the polar 
emission angle in the median plane ip, the planar trajectory length l xy from the SDC + CDC particle tracking, and 
the time of flight and z-coordinate (OH position) from the IH and OH scintillators: 

, 1 (,up 1 ,down\ , 
1 = 2 V'OH + t OH ) - HH 

Z=\{t d OH n -t a J H )Ve f f (6) 

The up-down indices correspond to the timing measurements at the two ends of the OH, and v e ff is the effective light 
transmission velocity in the scintillator material. These yield the primary observables (Figs. ki],p 



hi P = tan 



P = Pxy /cost 

I = l xy /cos6. 
/3 = l/ct 
m = p/Pjc 



z 



6 SC = cos 1 ( cosip I (7) 
V P J 

where B& ip is the out-of-plane dip angle, p and I the total momentum amplitude and three-dimensional trajectory 
length, and 9 SC the scattering angle with respect to the incident beam (Fig. ^). A left-right asymmetry noted in the 
scattering angle SC spectrum (Fig. ||d) has been reproduced in the MC simulations. 
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A coincidence of two charged particles, one on either side of the photon beam, signifies the occurence of an event 
(Section |[I C| ). A series of tests and cuts in the data set subsequently eliminate all but the ir + ir~ pairs. In particular, 
first the time-of- flight versus planar momentum spectra are used for particle identification (PID, Fig. ||). The great 
majority of events including a proton or are thus discarded. Cuts on the tagger TDC spectra reject events induced 
by spurious photons, not corresponding in timing with the beam pulse. Last, pairs with low-conhdence tracks (large 
RRMS), or whose vertex falls outside the target area (Fig. ^) are eliminated, thus completing a first-level selection 
based on directly measured observables. 

For the two-track events that have cleared the tests above, and have been identified as tt + tt~ , additional kinematical 
observables are calculated. At this stage, the few surviving two-track events involving a proton or e ± that had been 
previously misidentificd as tt + it~ by the PID cuts are eliminated as well. The calculated observables include the 
dipion invariant mass m„.+„.-, the laboratory production angle of the dipion system O/a/, the missing mass va miss 
and momentum p m i SS , the 7r + -7r _ laboratory opening angle i? 7r + 7r -, and the 7r + production angle in the dipion center 



of mass 0* + , employed as variables in the MC fitting procedure (Section |IV[ ) . 

Among these observables, the production angle for either one of the two pions in the dipion center-of-mass frame, 
for example 0* + , is singular as a direct experimental observable which, without the aid of simulations or assumptions, 
points to the presence of the p° production channel well below the nominal threshold energy. This is discussed next. 

C. The J = 1 signature of the p° 

Among the dominant production channels participating in tt + tt~ photoproduction in the E 7 ~ 1 GeV region 



(Section IVA). the p° — > tt + it~ channel alone results in the two pions being produced at a single reaction vertex 
with the characteristic J = 1 angular correlation from the decay of the p . In the dipion center-of-mass frame 
this translates into a pure cos 2 #* + distribution, where 6** + is the ir + production angle with respect to the dipion 
momentum, the direction defined by the latter in the laboratory frame. On the assumption of a slowly-varying 7r + 7r~ 
background interfering with the p° amplitude, the angular distributions are expected to be symmetric around 6* + 
= 90° for dipion cm energies near the mass of the p° meson, where the p° — » 7r + 7r~ amplitude peaks. Away from 
the p a mass, the resonant amplitude vanishes, and the background processes dictate the shape of the spectra [ j33| . 
Thus, above and below the /?°-meson mass, the angular distribution is expected to regain a quasi-isotropic shape due 
either to the uncorrelated pions produced at two or more reaction vertices, this being the case for the majority of 
the participating background processes, or from s-wave correlated pions, possibly produced from the decay of the a 
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meson (Section [IV A| ). 



This technique, of p '^ identification via the study of the pion-scattering angle in the dipion cm frame, has been 
extensively used in many previous analyses (e.g. Ref. |mJ). The 9* + distribution spectrum, based on the above, is 
expected to be well-described as A + B cos 2 #* + , in the vicinity of the dipion invariant mass matching that of the p° 
meson, where the p° — > 7r + 7r~ amplitude peaks. 

The 7r + 7r _ events in the range of 400-800 MeV/c 2 in dipion invariant mass have been divided in four 100 MeV/c 2 
bins, which is the finest binning allowed by the data statistics. An additional cut, determined from the MC simulations 
of the TAGX p° detection efficiency and kinematical considerations, eliminates those tt + tt~ events with too small an 
opening angle to have been the outcome of back-to-back production at a single reaction vertex (see Section [V Bj ). 



The latter cut results in a further 9-10% reduction in the total number of tt + tt~ events in the 400-800 MeV/c 2 dipion 
invariant-mass region, affecting only events from background processes, effectively boosting the p° — > ir + ir~ amplitude 
relative to the background (Fig. ||). 

The 500-600, and 600-700 MeV/c 2 regions (Figs. ^]b,c) clearly demonstrate the J = 1 fingerprint of the p° meson 
decay. The deviation from cos 2 0* + toward ±1 is reproduced in MC simulations of the p° — » tt + it~ process, and it 
is shown to be the effect of the TAGX detection efficiency, stemming from the two-track detection requirement (see 
Section IVB and Fig. [?]). The quasi-resonant p° — > 7r + 7r~ amplitude over the 500-700 MeV/c 2 dipion invariant-mass 



range points to a substantially reduced p° mass beyond the trivial apparent lowering, which is the artifact of probing 
the lower tail of the p° mass distribution with low-energy photons (Section |VB|) . 



IV. SIMULATIONS 

The MC simulations constitute an integral part of the data analysis by determining the process-dependent detec- 
tion efficiencies of the spectrometer, guiding the assignment of the weight to each of the contributing production 
mechanisms, and, ultimately, leading to the extraction of the medium-modified p° mass. 

The steps involved in the simulations and fitting algorithm can be outlined as follows: 



Eleven individual 7r + 7r production channels are coded into MC generators (Section [TV A ). These eleven pro- 



cesses are considered to account for the full ir + ir photoproduction yield in the 7 + 3 He reaction. Twelve 
distributions of six kinematical observables, with cuts aiming to separate the bound 3 He from the dissociated 



ppn final state, are simulated for each production channel and each of four AE 7 energy bins (Section IV 
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The analysis of the MC events is identical to that of the experimental data, and yields the process-dependent 
spectrometer acceptance (Section [V B| ). 



The simulated spectra for nine of the above elementary processess (Section IV A), including background and p° 



production channels with the nominal m p o=770 MeV/c 2 mass, are combined. The twelve spectra of each process 
are adjusted with a common strength parameter within each of the four AE 7 bins before being added together. 
Subsequently, all twelve simulated spectra are fitted simultaneously to the corresponding experimental ones, 
yielding the nine strength parameters independently for each AE 7 bin. From the latter, and the spectrometer 
acceptances, total cross sections are extracted for each of the nine production processes, and compared with 
independently established ones. Adjustments to the starting values and fitting constraints are made in iterative 
steps until satisfactory agreement is reached. 

The procedure is repeated for all eleven production channels, including the addition of two more processes 



(Section IV A) with a modified p mass m* in the range 500-725 MeV/c , but common for both the 3 He and 



breakup ppn final states. The p mass corresponding to the best fit for each AE 7 bin is quoted in this report 



as the medium- modified m* mass (Section VA) 



• Exploratory fits are attempted, decoupling the He and breakup ppn final states with respect to m* , as well 



as modifying the width r* (Section VB) 



The principal aspects of this algorithm are elaborated below. 



A. Production channels 



Several mechanisms are known to contribute to tt + it~ photoproduction. Recent experiments for photon energies 
between 450 - 800 MeV at MAMI, using the large-acceptance spectrometer DAPHNE and high-intensity tagged 
photon beams, and in the range 1 - 2.03 GeV with the SAPHIR detector at ELS A pq], have provided accurate 



measurements of the reaction jp — > tt + tt~p [|36 37 3S|]. These have motivated several theoretical models, which 
concur in their interpretation of the data as n + ir~ photoproduction predominantly through the 7rA(1232) — > tt + tt~ N, 
and the A^*(1520) — > ttA — > tt + tt~N channels |39|,[l0 41|. In the nuclear medium, the propagators of baryonic 



resonances require renormalization, and, in addition, many-body effects caused by pion rescattering (FS1) are known 
to interfere with the lowest-order reaction mechanism of two-pion photoproduction on the nucleon |42j] . These medium 
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modifications affect both the strength and the peak position of the cross-section spectra for the various interfering 
channels, relative to the corresponding processes on a free nucleon. Nonetheless, the A(1232) and N*(1520) resonances 
remain the leading channels in photon-induced reactions in the nuclear medium, as has recently been verified by total 
photoabsorption cross-section measurements on nuclei |l3|]. In the latter, substantial contributions were also attributed 
to the nucleonic excitations Pn(1440) and Sn(1535), primarily, which largely overlap with the N*(1520) resonance 
in medium-modified mass and width. 

The double- A is another channel that has been verified in photoabsorption measurements on the deuteron p4] , ^5[ , 
a process that has also been modelled theoretically ^(|. The photon is absorbed on two nuclei, exciting A(1232) 
resonances, in the reaction jNN — > AA, which subsequently decay to produce tt + tt~ pairs. 

In addition, thrcc-pion 7r + 7r~7r° production, associated with 3 He disintegration, is kinematically feasible in the 
energy range probed by the present experiment. However, the limited out-of-plane detector acceptance coupled 
with appropriate missing-mass cuts minimize the contributions of this mechanism. The experimental measurements 
available are sparse for this process in the energy regime of interest [[§6 47 4§| . 

Other possible contributions to the background 7r + 7r~ count, which, however, were not found to improve the quality 
of the fit and are presently not included in the simulations, may come from non-resonant three-, four-, and five-body 
phase space, corresponding to the 3 He remaining intact, or breaking-up into dp and ppn respectively. These multi-body 
phase-space processes are governed solely by energy and momentum conservation, each with a constant transition 
matrix element p9[ , and, loosely speaking, accomodate in an average sense all the remaining possible production 
channels which are too weak to be individually identified. 

The contributions of the mechanisms discussed so far have been previously considered in MC simulations, in 
connection with TAGX it + it~ photoproduction data |pl| , p3| , p5| . In Ref. |23| in particular, where it + it~ photoproduction 
via the p° channel was first considered, the background processes 

i) An(NN) sp 

ii) N*(1520)(NN) sp -> A(1232)7r(iVA% 
7 + ;i He -> { m) N*(l520)n(NN) sp 

iv) AAN sp 

v) ppnir + ir~ ir° 

were included in simulations of non-/?° tt + tt~ contributions, as well as final-state interactions (FSI) following the p° 
decay (see process ix) below). The index sp signifies spectator nucleons. The empirical values from Ref. j43| were used 



ppnir + TT 



1G 



for the A(1232) mass and width and for the N*(1520) mass, but the fit improved with the N*(1520) width doubled 
relative to Ref. p3| . This ad hoc increase effectively incorporates the near-by resonances Pn(1440) and Sn(1535), 
which largely overlap with the N*(1520), but cannot be resolved within the sensitivity of the data. Alternate fits were 
performed with the N*(1520) replaced by the Roper N*(1440) and including five-body phase space. The two methods 
yield comparable masses for the p , but the former is preferred as it results in a better fit. 

Additional improvements in the fitting procedure, relative to the analysis of Refs. pl"| , p3| , p5[ , include the modification 
of the MC generators to account for the angular momentum of all p a — > 7r + 7r~ and intermediate A-resonance channels. 
Furthermore, motivated by recent tttt phase-shift analyses which increasingly show evidence of s-wave contributions 
from the a meson, the quasifree cr-decay channel 



vi) 7 + 3 He — > appn — > tt + tt 3 He 



has been added, with the a mass and width parameters from Ref. [f50|. 

Last, p° — > 7r + 7r~ photoproduction has been simulated by means of five distinct generators, namely 



7+ 3 He 



vii) p° + 3 He — > 7T+7T 3 He 



viii) p°ppn — > 7T + 7r ppn 



ix) p°ppn — > 7r + 7r ppn + (FSI) 

x) p° + 3 He -> tt+tt- 3 He 



m p o = 770 MeV/c 2 



xi) p°ppn — » 7T + 7r ppn 



500 MeV/c 2 < m* < 725 MeV/c 



where to channels vii)-ix) and x) — xi) are ascribed p° decay outside, and inside the nuclear medium, respectively, 
the latter probing the medium effect on the p° mass. The breakup channels have the reaction taking place on a single 
nucleon, subject to its Fermi motion, with the remaining two nucleons as spectators. Final state irN interaction (FSI) 
with one of the two spectator nucleons is included in channel ix). 
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B. Fitting procedure 



In modelling tt + tt photoproduction on nuclei, the distributions of five kinematical observables (Section [II B) were 
simultaneously fitted to the data in Refs. HJ23,E3]. These are 



1. the dipion invariant mass m n + w - 

2. the laboratory production angle of the dipion system 0jm 

3. the missing mass m m i SS 

4. the missing momentum p m iss 

5. the 7r + — 7T _ laboratory opening angle 



to which one additional kinematical observable has presently been added (Section III C ) , namely 



6. the 7r + production angle in the dipion rest frame cos9* + . 

Moreover, in Ref. J2lj it was determined that dividing the data sample in AE 1 — 80 MeV bins provided the 
optimal compromise between the presumed constancy of the reaction matrix elements, implicit in the MC simulations 
which only depend on the kinematics, and the requirement of sufficient statistics. The AE 7 partitioning of the data is 
necessary in order to account for the varying energy dependence of the tt + it~ cross sections from each of the individual 
production mechanisms. The 80 MeV binning in E 7 has been retained, resulting in four sectors of the data sample 
from 800 to 1120 MeV, to be referred by their respective central E 7 values (840, 920, 1000, and 1080 ±40 MeV). 

The addition of the cos#* + distribution, though not noticeably affecting the overall quality of the fit, nonetheless 
provides an additional physical constraint which aids the MC fitting algorithm to converge to a more realistic solution. 
In particular, this kinematical observable uniquely captures a characteristic feature of the contributing mechanisms, 
which may be classified into three types according to their respective dependence on cos#* + (Fig. |^): 

• The channel of interest, diffractive p° — * 7r + 7r~, is unique in producing two p-wave correlated pions. The 
spectrometer response to this mechanism is consistent with the J — 1 dependence, and the deviation from 
the anticipated cos 2 0* + distribution towards ~ ±1 reflects the acceptance cut, stemming from the kinematical 
conditions, set-up geometry, and the two-pion detection requirement (compare the solid curve of Fig. |7] with 
Figs.|b,c). 



The background hadronic channels i) — iv) of Section IV A produce two uncorrelated pions at two or more 



reaction vertices. The angular correlations of these pions are averaged out over Air sr in simulations, resulting 
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in featureless cos#* + spectra. The spectrometer geometry, however, suppresses the pion acceptance away from 
cos#* + =0 (see e.g. the dashed curve of Fig. ^| for the single- A production channel). This is a consequence of 
the fact that channels i) — iv) involve the decay of intermediate baryonic resonances, accompanied by energetic 
nuclcons. 



Three-pion production and the quasi-elastic a process, v) and vi) of Section IV A| , are characterized by featureless 



cos6** + acceptances, as no energetic nucleons are emitted (dotted curve of Fig. ^). 

The combination of improvements relative to the analysis of Ref. |23j | , namely, accounting for the angular momentum 
in the A and p° channels, and imposing additional physical constraints via the new kinematical observable cos#* + , 
resulted in a more accurate treatment of the process-dependent spectrometer acceptances. 

The data have been subjected to two additional cuts, one of which enhances the p° relative to the background 
channels, and the other facilitates the separation of the bound 3 He from the breakup ppn final states. The former is a 
7r + 7r~ opening-angle cut determined from MC simulations, namely 70° < < 180°, eliminating two-pion events 

that could not have been produced back-to-back from the p° decay (Fig. ||a). This cut is most effective at the higher 
end of photon energies covered by the experiment, where the p° identification becomes difficult by an increasingly 
deteriorating yO°-to-background ratio with increasing E 7 . The latter cut (Fig. ||b) separates events with missing mass 
in the proximity of the target mass mi He ~ 2.8 GeV/c 2 (i.e. 2700 MeV/c 2 < m m i SS < 2865 MeV/c 2 ), from those 
associated with the breakup of the target nucleus to ppn in the final state (2865 MeV/c 2 < m m i SS < 3050 MeV/c 2 ). 
The combination of the two types of cuts is applied to three of the kinematical observables, resulting in six additional 
spectra, besides the unselected tt + tt~ distributions 1-6 enumerated earlier. These are: 



7. rriK+ir- 

8- Pmiss 

9. cos0* + 



10. m v +T t - 

11- Pmiss 

12. cos9*. 



> p° — enhanced low missing — mass data 



> p° — enhanced high missing — mass data 



While the contributing p° photoproduction mechanisms may not be experimentally distinguishable (Section ||), the 
cuts aim at partitioning the data in biased samples favoring processes which are more prone to probing either the 
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nuclear mean-field effect (i.e. coherent and exclusive quasifree photoproduction via the distributions 7-9), or a possible 
nucleonic effect (i.e. the breakup quasifree channel via the spectra 10-12). The three kinematical observables which 
were subjected to the cuts, namely, the dipion invariant mass, missing momentum, and dipion-cm ir + production- 
angle, were selected empirically from the kinematical observables 1-6 as more sensitive to the p° mass, and therefore 
more susceptible to possible medium effects. 

With a range of m* values, traversing the region 500-725 MeV/c 2 , and each value kept common for the unselected, 
as well as the 3 He and ppn selected data, the twelve simulated spectra were fitted simultaneously to the corresponding 
experimental ones, by minimizing a standard x 2 function with the strengths of the eleven individual processes i) — xi) 
(Section |TV A ) as the fit parameters. The four AE 7 = 80 MeV bins were fitted independently. The optimal m* for 



each bin is that corresponding to the minimum value of the \ 2 function (Section V A). 
In summary, the MC fitting procedure satisfies the following requirements: 

• The AE 7 = 80 MeV binning restricts the energy dependence of the participating processes to the narrowest bin 
possible without loss of sufficient statistics. 

• The six kinematical observables utilized in the fitting are complementary, and account for different physical 
attributes of the data sample. This imposes far more stringent constraints than an analysis based on only the 
invariant mass distribution, as in the case of the CERES data || . 

• The simultaneous fitting of selected and unselected data aims at isolating a strong signal from data samples 
most responsive to possible p° mass modifications, while also incorporating in the fit the bulk of tt + tt~ events 
produced in processes less sensitive to such effects. 

The outcome of this procedure is the medium-modified p° meson mass. 



V. RESULTS 

A. Quality of fit and uncertainties 

Beyond the statistical and other experimental uncertainties which are folded into the calibration and analysis of the 
measured tt + tt~ yields, additional sources of uncertainty are generated by the MC simulations and fitting algorithm, i) 
The MC event generators depend exclusively on kinematical parameters, neglecting other aspects of the interaction. 
As an example, the channel 7 + 3 He — > (np) sp ATr is modelled as 7 + 3 He — > (np) sp A ++ 7r~, normalized by an 
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isospin scaling factor to account for the remaining hadronic charge states. This introduces an uncertainty in the 
amplitude of this process, similar to which are also present in other background hadronic channels, ii) Independent 
total cross-section measurements of the constituent reactions, serving to anchor their relative strengths in the full 
7r + 7r~ production process, are sparse (e.g. the quasi-elastic a channel). The strength for some of the individual 
channels was inferred from indirect sources, iii) Additional quasifree channels and independent measurements to fix 
their strength are required in order to extract precision cross sections for the background hadronic channels from the 
data, processes which merit attention on their own behalf. This may become possible following the analysis of three 
charged-particlc events from TAGX experiments, currently in progress. 

Despite the caveats above, the medium-modified p° masses extracted from the MC fitting procedure have remained 
remarkably stable with respect to variations of the strength parameters of the constituent reactions, within each of 
the four AE 7 bins. This is all the more significant considering that the data are fit simultaneously for twelve spectra 
of six kinematical observables. In conjunction with the direct J — 1 fingerprint discussed earlier, the insensitivity 
of the fit, within physical constraints, adds confidence to the premise that the medium-modified p°-meson masses 
extracted from the MC simulations reflect genuine features of the data sample. 



Following the procedure discussed in Section IV, several MC fits have been performed with m* taking on values 
from a mesh in the range 500-725 MeV/c 2 (Fig. ^J). The steepness of the x 2 - v s-m* curve is indicative of the sensitivity 
of the data sample to p° mass modifications, within each of the four Ai? 7 bins. Using this as a qualitative criterion, 
the fitting of the 840 and 920 MeV samples, and to a lesser extent of the 1000 MeV sample, have converged to a 
"best" m* , whereas the fit for the 1080 MeV bin is essentially insensitive to variations of m* (Fig. ||). 

For each of the 840, 920, and 1000 MeV bins individually, the MC fits (dark circles, triangles and squares respectively 
in Fig. ||) yield the best m* and an estimate of its uncertainty. This is achieved as follows: a) The uncertainty a x 2 is 
assumed common within each bin. This is justified by the fact that both the data set and the reaction matrix elements 
used by the MC algorithm are common within each bin. b) A polynomial expansion of the MC x 2 about (m* -mo), 
with mo in the proximity of the apparent minimum, is assumed to describe well the dependence of the former on the 
latter. Subsequently a new \' 2 minimization yields the rank and the coefficients of the polynomial. A third order 
polynomial gives the best result for the two lower-energy bins, and it is also assumed to provide the best description 
of the parent population for the 1000 MeV data sample. This procedure also leads to conservative estimates for a x 2 , 
in particular 0.019, 0.013 and 0.017 for the 840, 920, and 1000 MeV bins, respectively, c) The polynomial coefficients 
and a x 2 estimates are used to derive the optimal m* and an estimate of its uncertainty for each bin. These are: 
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• m* = 642±40 MeV/c 2 , 800 MeV<E 7 <880 MeV 

• m* = 669±32 MeV/c 2 , 880 MeV<E 7 <960 MeV 

• m* = 682±56 MeV/c 2 , 960 MeV<E 7 <1040 MeV 

The improvement in the medium-modified over the unmodified p° mass fits is evidenced to varied degrees in all the 
fitted spectra. The order of the improvement is illustrated for the full (unselected) data in the dipion invariant-mass, 
7r + -7r~ laboratory opening-angle, and missing-momentum distributions for the 840±40 MeV and 920±40 MeV bins, 
which are most affected by modifications of the p° mass (Figs. |Ic|Jl"l"| ). 

The dipion invariant mass spectra are not the most sensitive to variations of the p° mass, among the twelve 
distributions of six kinematical observables that were employed in the fit. This is seen, for example, in comparing 
the improvement between the unmodified- and modified-mass dipion invariant-mass with the opening-angle and 
missing-momentum spectra (Figs. |l^Jll|), where the latter two observables are seen to display a greater sensitivity. 
This underlines the advantage of a fitting procedure that utilizes multiple complementary kinematical observables as 
opposed to only the invariant mass, thus capturing additional attributes of the physical process, and resulting in a 
more realistic analysis. 



B. Discussion 



The J—l angular momentum signal, discussed in Section III C (Fig ph, as well as the dipion invariant mass spectra 



from the fit (e.g. Figs. ^a,d for E 7 =840 MeV), would appear to indicate medium-modified p masses which are 
actually lower than the values quoted earlier (e.g. compare m* o =642±40 MeV/c 2 for the 840 MeV bin with the 



centroid suggested by the dashed curve of Fig. |lO|a). This apparent discrepancy is misleading, and has its origin in 
a trivial effective mass lowering driven by phase-space, most pronounced at lower photon energies. This is due to 
the fact that, at low photon energies, only the lower wing of the p° mass distribution is kinematically accessible (e.g. 
Fig. [l^ for E 7 =840 MeV and the p°ppn final state). The shape and the centroid of the mass distribution for the 
resulting p° mesons is primarily dictated by the kinematics, and to a lesser extent by the spectrometer acceptance 
for the particular p°-producing process. This effect is implicit in the MC generated p° spectra, and it is by no means 
sufficient to account for a mass lowering of the order explicitly observed in the experiment. This is manifest, for 
example, in comparing the lower solid curves of Figs. [To|a,d - for which p° production is dominated by the ppn final 
state - with the J=l fingerprint of the p° in Fig. |^. The loci of the former curves, consistent with the mass distribution 
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indicated by the histogram of Fig. [l2| a for E 7 =840 MeV and the nominal p° mass, are far too high to be compatible 
with Fig. ^, with a p a signal peaking in the range of 500-600 MeV/c 2 . In contrast, the histogram of Fig. |l2] b for a 
lowered p° mass, and the dashed curve in Fig [To|a, are consistent. 

Alternate fits were also attempted, to possibly discern additional features from the 840 and 920 MeV bin samples. 
Decoupling the exclusive 3 He from the breakup ppn final states with respect to m* yielded identical masses for 
the 920 MeV bin, and a somewhat smaller mass for the latter relative to the former channel for the 840 MeV bin. 
The resulting improvement in the quality of fit, however, is within the uncertainty estimate of o~ x 2. Fits were also 
performed with r* fixed to the predicted width at half nuclear density extracted for 3 He from Ref. Jl3|], about double 
the free width T p a . A sizeable improvement in \ 2 was noted with the modified width for any mass, compared with 
the free mass and width case. In all cases, however, the \ 2 is 5-10% larger with r* than with r p o. Moreover, the 
preferred m* is no higher with the modified than with the free width. In summary, these exploratory fits verify 
the preference in a reduced p° mass, but are inconclusive, within the sensitivity of the data, as to whether a width 
modification is in addition supported. 

The absence of a conclusive p° mass-modification dependence from the 1000 and 1080 MeV bins, in contrast to 
the strong signal from the 840 and 920 MeV bins, is telling as well. Whereas the former are more prone to probing 
the exclusive channel, and therefore the nuclear mean-field medium effect, the latter are deeper into the subthreshold 
region, dominated by the breakup channel, and probe distances shorter than the nucleonic radius. To illustrate the 
range of the p° processes in different energy regions, we consider the mean decay length Iq of the p° 

'o = Pto = — — 
rripO l p oc 



m p o7c 

h 

to = 7^0 = , (8) 

1 p 

with the nominal mass m p o and width r p o for the p° meson, in the rest frame of either the interacting nucleon, or 
the 3 He nucleus, for the breakup or bound final state (Fig. |l3|). At low photon energies (e.g. E 7 =840 MeV, Fig. |l3|a) 
the ppn channel dominates p° photoproduction. This is verified by the fact that the low missing-mass selected data 
represent less than 10% of the total unselected events, whereas the high missing-mass selected data contribute the great 
majority of the p a events and their corresponding distributions are generally consistent with the full (unselected) data. 
With the p° mesons produced on the nucleon, and a mean decay length of under 0.5 fm for a substantial portion of 
them, there is a large overlap of the volume traversed before their decay and the nucleonic volume. To the extent that 
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the medium modifications of the vector-meson properties depend on the density of the surrounding nuclear matter, 
it is conceivable that the induced medium effect be dominated by the large densities encountered in the interior of 
the nucleoli, and consequently be far more pronounced than a medium effect inducted by nuclear fringe densities. 
At higher photon energies (e.g. E 7 =l GeV, Fig. |l3|b), both the ppn and 3 He final states contribute and the mean 
decay length increases. Whereas nucleonic densities become increasingly inaccessible, a large overlap of the volume 
traversed by the p° before its decay with the 3 He nucleus may still induce a weaker nuclear medium effect. Moreover, 
the possibility of a medium-induced increase in the width r p o is in favor of shorter decay lengths (Eq. |^), and therefore 
further increases the likelihood of accessing regions of large densities either in the nucleon, or the nuclear core. The 
present results may therefore be suggestive of a moderate medium effect in the realm of the nuclear mean field at 
near-threshold photon energies, probing normal nuclear densities, turning to a drastic reduction of the p° mass in the 
deep subthreshold region, where the p a decay may be induced in the proximity of the nucleon. These implications, 
however, need to be further investigated and verified by higher precision and large solid angle experiments p§| ]. 



VI. CONCLUSIONS 

In summary, the kinematics and final state of the 3 He(j,ir + ir~)X reaction have been studied with the TAGX 
spectrometer and a tagged photon beam in the subthreshold p° photoproduction region. The bound 3 He and breakup 
ppn components of the p° channel have been investigated, aiming at the distinction between a nuclear and a possible 
nucleonic medium effect on the p° mass. The p° channel has been aided by the inherent selectivity of the TAGX 
spectrometer to coplanar ir + ir~ events, and by the choice of 3 He as a target with minimal FSI effects without 
suppression of the p° amplitude. In any case, the size of target has little effect in the subthreshold region. The 
J = 1 fingerprint of the p° has been observed in the dipion invariant mass region 500-700 MeV/c 2 , pointing to 
a substantial reduction beyond a trivial phase-space governed apparent lowering. This has been verified by MC 
simulations, incorporating the exclusive and breakup p° channels, the latter both with and without FSI, as well as 
background hadronic channels and s-wave correlations. The extracted p° medium-modified masses, 642±40 MeV/c 2 , 
669±32 MeV/c 2 , and 682±56 MeV/c 2 for the 840, 920, and 1000 MeV data bins, suggest a strong medium effect in the 
deep subthreshold region, requiring large densities that are incompatible with a nucleus as light as 3 He, but that are 
more consistent with the probing of the nucleonic volume. The pattern of decreasing sensitivity with increasing photon 
energy, with the 1080 MeV bin showing no evidence of a p° mass-modification signal, hints a moderate mean-field 
nuclear effect at near-threshold energies. The simulations are inconclusive regarding a medium-modified p° width. 
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Further analysis currently in progress for photoproduction on heavier targets ( 4 He, 12 C) and for background con- 
tributions (e.g. the An channel) from TAGX experiments, as well as future experiments planned for TJNAF p8[ , 
may more accurately isolate the p° — > ir + ir~ channel from the background processes, and shed light on the nature of 
the medium modifications on light vector- meson properties at the interface of hadronic and quark matter. 
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FIGURES 




FIG. 1. Components of the TAGX spectrometer: Radially outwards from the target vessel at the center are seen the IH, 
SDC, CDC, and OH. The TAGX coordinate frame has the ai-axis pointing in the direction of the photon beam, and the 5-axis 
pointing outward from the page in the direction of the 0.5 T magnetic field of the dipole analyzer magnet. Veto counters (not 
shown in the figure) are positioned along the OH arms in the xy-plane, for e + e~-pair rejection, mainly at forward angles. 
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FIG. 2. The vertex position is plotted, reconstructed from the 8-layer CDC data only (left panel), and from the 9-layer 
CDC and SDC data (right panel), for a set of empty-target runs [Q. The evident improvement in the resolution is discussed 
in the text. The components of the target vessel, visible in the figure, have been discussed elsewhere [ p9|j3C| ]. The inner ring 
corresponds to the Mylar wall containing the liquid 3 He target, a cylinder 50 mm in diameter. A circle indicating the position 
of the target container has been drawn, centered at (x,y)=(2.8 mm,-1.5 mm) in the TAGX coordinate frame. 
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FIG. 3. The x-coordinate of the vertex position after background subtraction is shown, along with the target-vessel walls 
indicated as lines. Only ir + ir~ events from this target region are considered in the analysis. 
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FIG. 4. The TOF-vs-planar momentum spectra, displaying proton (upper- most), ^ (in the box cuts), and e (adjacent to 
the p xy — axis) bands, are used for particle identification. 
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FIG. 5. Some of the kinematical observables for tt + tt~ events, described in the text, are shown. The sign of the momentum 
coresponds to the pionic charge. The sign of the scattering angle 8 S c depends on whether the track registered to the left (>0) 
or right (<0) of the beam (see Fig. [l]). 
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FIG. 6. The cos#* + distribution captures the J = 1 signature of the p° — > tt + tt~ decay. Panels b) and c) display the 
A+Bcos 2 £T + dependence (dashed curve), expected on the basis of the J = 1 angular momentum, and the deviation towards 
±1 is due to the spectrometer acceptance. 
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FIG. 7. Simulated events from different processes, taking into account the spectrometer acceptance, have one of three 
characteristic cos#* + profiles. The simulations shown are differential cross sections normalized to unity at E 7 =920 MeV, and 
they display the cos#* + dependence for three representative reactions: a) 7 + [i He — > p°ppn — > -K + ir~ppn with the nominal m p o 
(solid curve), b) 7 + 3 He — ► 7r~A ++ pn — > -K + -K~ppn (dashes) and c) 7 + 3 ffe — » appn — ► 7r + 7r _ ppn (dots). 
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FIG. 8. The effects of the a) tt + -tt~ opening-angle and b) missing-mass cuts, indicated as arrows in the figure, are shown 
for the respective simulated spectra of the p° 3 He, p°ppn and 7r~ A ++ pn channels (solid, dotted and dashed lines respectively). 
The calculations are for E 7 =920 MeV and the integral of each distribution is normalized to unity. 
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FIG. 9. The dependence of the \ 2 function of the MC fits on the variation of the p° mass is shown for the four AEj — 
MeV bins. The MC calculations are indicated as points, and fits to third-order polynomials as curves. The fitting procedu 
yields the best m* per A£ 7 bin, and estimates of the respective uncertainties. 
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FIG. 10. The dipion invariant-mass (panels a, d), -k + -tt~ laboratory opening-angle (panels b, e), and missing-momentum 
(panels c, f) spectra are shown for the full (unselected) data at E 7 =840 MeV. The solid squares are the data, showing only one 
representative error bar for clarity. The remaining error bars are comparable and have been included in the fitting. The MC fits 
shown are: full reaction (upper solid lines), non-p° background processes (dotted lines), p° contributions with unshifted mass 
(lower solid lines), and p° contributions with m* =650 MeV/c 2 (dashed lines). The latter are included only in the calculations 
represented by the top panels (a, b, c). The improvement from bottom to top reflects the effect of including the p° medium 
modifications. This improvement is more pronounced for the p° selected data set, not shown. 
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FIG. 11. The data and curves shown are as in Fig. [u], but for the E 7 =920 MeV bin, and m* =675 MeV/c 2 . It is stressed 
that the fitting was done independently for each of the four AE 7 =80 MeV/c 2 bins, as discussed in the text. 
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FIG. 12. The trivial p° mass shift due to kinematic phasespace limitations is illustrated for the quasifree 
7 + 3 He — ► p°ppn — > TT + Tr~ppn process at E 7 =840 MeV, a) for the nominal and b) for a lowered p° mass. The MC gen- 
erator randomly selects p° masses from the Gaussian distributions (solid curves). The histograms (dashed curves) are the 
resulting m p o spectra whose integrals are normalized to unity, after verification that the reaction is kinematically feasible with 
the selected mass and that the photoproduced pions are accepted by the spectrometer. 
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FIG. 13. The mean decay length lo distribution, of the p° — > -k + tt~ decay, is indicative of the distances probed for vec- 
tor-meson modifications in this experiment. The lo spectra are illustrated for a) E 7 =840 MeV and b) E 7 = 1000 MeV. The 
histograms are for MC generated events accepted by the spectrometer in the 7 + 3 He — > p°p{pn) ap — > tt + tt~ p(pn) sp (solid 
curves) and 7 + 3 He — > p° 3 He — ► tt + tv~ 3 He (dashed curves) reactions, with their integrals normalized to unity. The mean de- 
cay length is for the former reaction in the rest frame of the participating proton (the two remaining nucleons being spectators), 
and for the latter in the rest frame of the 3 He nucleus. 
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